Exact density-functional potentials for time-dependent quasiparticles 
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We calculate the exact Kohn-Sham potential that describes, within time-dependent density- 
functional theory, the propagation of an electron quasiparticle wavepacket of non-zero crystal mo- 
mentum added to a ground-state model semiconductor. The potential is observed to have a highly 
nonlocal functional dependence on the charge density, in both space and time, giving rise to features 
entirely lacking in local or adiabatic approximations. The dependence of the non-equilibrium part 
of the Kohn-Sham electric field on the local current and charge density is identified as a key element 
of the correct Kohn-Sham functional. 
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Density functional theory (DFT) in the local density 
approximation (LDA) has proven highly effective in the 
ground-state calculation of electronic structures and en- 
ergies of a diverse range of systems [T|. It was natu- 
ral therefore to extend the approximation to the adi- 
abatic LDA (ALDA) within a time-dependent density- 
functional-theory (TDDFT) description of excited mat- 
ter. However, while the ALDA is effective in certain 
classes of excited systems [2], it is shown to break down 
in a large number of physical scenarios of universal inter- 
est, for instance the absorption spectra of semiconductors 
and insulators [3], the polarization of systems subject to 
an external electric field [H [5] , and the prediction of cor- 
relation energies [3] and double excitation spectra [7]. 

One important class of non-equilibrium problems for 
which accurate TDDFT functionals are needed is the 
quantum transport of charge through matter. Histor- 
ically, quantum treatments of electronic transport re- 
duce the current response for an applied potential across 
a device to a transmission problem, and employ non- 
interacting electrons. More accurate calculations make 
use of effective potential fields including Hartree, ex- 
change and approximate exchange-correlation terms. A 
natural step, then, is a Kohn-Sham (KS) TDDFT ap- 
proach in which the exact potential is a functional of the 
initial state and the time-dependent density. The need 
to better understand Kohn-Sham functionals, including 
those depending explicitly on the current density, beyond 
the adiabatic and exchange-only regimes in order to de- 
scribe quantum transport has been previously noted [5]. 
Maitra et al. [5] have also highlighted the importance 
of including the initial-state dependence, absent in adia- 
batic approximations, of the time-dependent Kohn-Sham 
potential in the desired functional. 

To advise the construction of appropriate function- 
als, we believe it is important to study the exact Kohn- 
Sham potential itself. The reverse-engineering of ground- 
state KS potentials via direct numerical optimization 
of characteristic parameters has provided insight into 
the nature of nonlocal functional dependence of, for in- 



stance, polarized materials ^ [5], while in the time- 
dependent regime there have been calculations of exact 
time-dependent exchange-correlation potentials and ker- 
nels for non-equilibrium systems such as doubly-excited 
helium |10j and small Hubbard chains subject to time- 
dependent electric fields 

In this paper we introduce a self-consistent algorithm 
for reverse-engineering the exact potential from known 
time-dependent charge and current densities. To this 
end, we consider systems which may be solved exactly by 
other means, and yet which contain important physical 
aspects of quantum transport. Our model system is an 
electron propagating through a infinite, one-dimensional 
semiconducting wire. The addition of a single electron 
to a ground-state device - in this instance the semicon- 
ductor wire - is the domain of quasiparticle theory, in 
which the wavefunction and energies are given by the 
Schrodinger-like equation: 



Vext + V H - E 
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(1) 

where v ex t is the external potential, vh the Hartree po- 
tential, E the quasiparticle energy and £ the nonlocal, 
energy-dependent self-energy operator. 

To focus on the physics arising from the nonlocality of 
the self-energy operator, we take the net local potential 
of the ground-state V ex t + Vh to be zero and employ 
a model operator that is nonlocal but real and energy- 
independent [41IT2"]: 



E(x,x') = 



f(x) + f(x') 



i\x-x'\) 



(2) 



where g(x) = exp (— (x/u>) 2 ) / ^/nw introduces the nonlo- 
cality of the self-energy, and f(x) — —F [1 — cos(27ra;/a)] 
is the local, periodic potential to which £ would reduce 
if w were small. We choose our operator to approximate 
that of silicon, with a self-energy amplitude of i*o = 4.1 
eV, lattice spacing a — 4 a.u. and nonlocal range w = 2 
a.u. We model the infinite wire as a chain of "atoms" each 
contributing one spinless electron to the ground-state. 
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The Kohn-Sham ground state is then described by a 
scalar potential which reproduces the model ground-state 
density exactly, which, since it is periodic, may be found 
by direct numerical optimization of its Fourier compo- 
nents. The resultant KS band structure (Fig. [T]) dif- 
fers from that of the nonlocal model, with a band gap 
that falls short of the quasiparticle gap by 0.31 eV [I]. 
The difference in band structure gives rise to a differ- 
ence in group velocity (v g = dE/dk) between the Kohn- 
Sham and quasiparticle eigenstates, which would prop- 
agate a wavepacket incorrectly and thus necessitates a 
time-dependent component in the KS potential. 

In contrast to typical approaches to nonequilibrium 
systems, we will prepare the system "in flight", taking 
an arbitrary instant of the quasiparticle propagation as 
our t — 0. The addition of a pure quasiparticle with a 
single Bloch wavevector to a ground-state semiconduct- 
ing wire places the device in the steady-state regime; 
in this paper we focus on an explicitly time-dependent 
system: the addition of a quasiparticle wavepacket with 
a non-zero crystal momentum that is sufficiently well- 
defined that the quasiparticle approximation can be ap- 
plied. Thus we construct the quasiparticle wavepacket 
as a weighted sum over a narrow range of solutions to 
the quasiparticle equation with a momentum expectation 
value (pqp) = 0.65 x 2n/a. 

Because of the initial-state-dependence of the exact 
time-dependent potential, it is common to prepare the 
Kohn-Sham system in its ground state, where the exact 
potential is a unique functional of the density, and, at 
some arbitrary time, induce an excitation with a chosen 
time-dependent potential or else excite it toward some 
desired state or density via optimal control |13| . How- 
ever, the KS wavefunction at some intermediate time be- 
tween its preparation in the ground state and a later ex- 
cited state of interest implies a prehistory without further 
explicit functional dependence on the time-dependent 
density that preceded it [S]. Thus we may choose an 
initial KS state that is in-flight [H] , as for the model sys- 
tem, first with the condition that it yields the same initial 
charge density, and then by using the transformation 

i/)(x,Q) «- e l6{x ^i;(x,0) (3) 

and optimizing the real field 9{x) to yield the correct 
current density. 

The quasiparticle wavepacket, upon being introduced 
into the supercell, propagates without the necessity of a 
potential difference across the wire. (The Crank-Nicolson 
method [TS] is employed to solve the time-dependent 
Kohn-Sham and quasiparticle equations explicitly.) Fig- 
ure [l] shows snapshots of the time evolution of the quasi- 
particle charge density riQp(x,t) = n(x,t) — nos{x) and 
current density jqp(x, t) as it propagates a short distance 
along the wire. For reasons that will become clear below, 
we restrict ourselves in this paper to small displacements 
of the wavepacket. 
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FIG. 1: (Color online) The charge density additional to the 
ground state (top) and the current density (bottom) of both 
the model and Kohn-Sham systems at 1 time-step (solid red) , 
75 time-steps (dashed green) and 150 time-steps (dotted blue) 
of At = 4 x 10~ 3 a.u. The quasiparticle (solid black) and 
Kohn-Sham (dashed blue) band structures are inset, with the 
gray shaded region marking the range of pure quasiparticle 
Bloch states included in the electron wavepacket. 



The algorithm outlined in Figure [2] is then invoked at 
each time-step to calculate the exact Kohn-Sham po- 
tential. At the heart of the algorithm lies the obser- 
vation that, if the instantaneous current density for a 
Kohn-Sham system in the adiabatic limit is incorrect 
by an amount Aj — j(x,t) — jKs(x,t), one can alter 
the effective vector potential AKs{x,t) to yield a better 
current density j KS (x,t) = j [^Ks] + A KS (x,t)n(x,t) 
(where jo[^f] is the standard single-particle current den- 
sity) without affecting the instantaneous charge density 
n Ks{ x it), since the adiabatic vector potential performs 
a transformation of the Kohn-Sham system of the form 
of Eq. 3. The algorithm demands that the time-step em- 
ployed in the time-evolution of the Kohn-Sham system 
must be small enough such that, if the correct KS cur- 
rent density is obtained for consecutive time-steps, the 
change in charge density is correct by necessity of the 
continuity equation dn/ partialt + V • j = 0. Since the 
history of the KS density must match that of the quasi- 
particle system, the same time-step (here 4 x 10~ 3 a.u.) 
is used for both. 

The vector potential, while useful procedurally (Fig. 
[2]), is known to generally have a nonlocal (in time) 
functional dependence on the charge density [IBJ. Of 
more universal relevance is the Kohn-Sham electric field 
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FIG. 2: The reverse-engineering algorithm. A self- 
consistent calculation of the exact (in the limit of At — > 0) 
time-dependent Kohn-Sham potential which is repeated for 
every time-step in the quasiparticle system history. The algo- 
rithm makes incremental corrections to the Kohn-Sham cur- 
rent density, with the correct charge density arising due to the 
continuity equation. Here, vks is the Kohn-Sham potential, 
$ is the Kohn-Sham Slater determinant, Aj is the current 
density error, and Ut,t-At is the unitary evolution operator 
which time-evolves from time t — At to time t. 



which may be defined as $Ks{x,t) — —9tAKs{x,t) — 
d x VKs{x-,t) which may then be implemented in the 
Kohn-Sham equations via Kohn-Sham vector and/or 
scalar potentials as appropriate. 

Snapshots of the exact additional Kohn-Sham scalar 
potential for the quasiparticle wavepacket system are 
shown in Figure [3] We observe first that the exact poten- 
tial is not periodic in space. For this reason, and noting 
that the quasiparticle is itself finite in space, only a fi- 
nite section of the Kohn-Sham infinite wire was modelled, 
specifically a chain of 20 "atoms" giving a total electron 
number of 21. We observe that, at all times, the charge 
densities on each side of the wavepacket region are that 
of the ground-state wire, while the potentials differ by a 
time-dependent constant. This ultranonlocal nature of 
the potential difference negates the possibility of the ex- 
act potential being approximated by a local functional 
of the charge, or indeed current, density. Such poten- 
tial steps in the exact time-dependent Kohn-Sham po- 
tential have been observed before in electron ionization 
calculations [17] and have been linked to the exchange- 
correlation derivative discontinuity |18j as the electron 
number of the system passes through an integer. 

Our results differ from previous calculations, however, 
in three respects. First, although the electron number 
in a finite section of our system will gradually increase 
from N to N + l and back again as the wavepacket passes 
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FIG. 3: (Color online) The additional (relative to the ground 
state) time-dependent Kohn-Sham potential Avks (top) and 
electric field ASks = —dvxs/d x (bottom) after 1 time-step 
(solid red), 75 time-steps (dashed green) and 150 time-steps 
(dotted blue) of At = 4 x 10 _3 a.u. The gradient is local- 
ized and more slowly-varying, making it more amenable to 
approximation than the potential itself. 



through it, the electron number of the system as a whole 
remains constant. Second, the potential step of |17| takes 
the form of a plateau which expands from the central re- 
gion as the ionized electron escapes the system, while, in 
our system, the ends of the wire are immediately placed 
in a potential difference due to the localised electric field 
of the added electron. Third, this potential difference is 
time-dependent, thus a section of the wire the wavepacket 
has already passed through will continue to react to elec- 
tron motion in principle infinitely far away. 

By comparing Figures [T] and [3j which feature snapshots 
at the same time intervals, one observes that this nonlocal 
aspect of the potential time-varies much more rapidly 
than the wavepacket charge or current density. Semi- 
classically, one can understand the necessity that such 
an electric field be time-varying: an electron subject to a 
constant electric field will acquire a constant acceleration. 
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Thus a static potential difference cannot reproduce the 
steady progression of the actual quasiparticle, and a time- 
varying field must be present. 

The additional KS potential (Figure [3] also exhibits a 
small barrier, localized with the quasiparticle, and mod- 
ulated by a periodic function with the periodicity of the 
crystal lattice. The barrier acts as a source of resistance 
to the current density for the Kohn-Sham electrons while 
the periodic part of the additional potential acts to 'tune' 
the local band structure, and hence the group velocity, 
of the Kohn-Sham electrons and, in conjunction with the 
rest of the barrier, ensures that the shape of the charge 
density evolves correctly. 

It is the combination of these three components of the 
exact time-dependent potential that provides the mech- 
anism for yielding the correct current density. In con- 
trast to the quasiparticle description, in which all of the 
current is carried by the added quasiparticle, a small 
but significant proportion (about 5%) of the Kohn-Sham 
current density is carried by the valence-band electrons 
which are subject to the same additional potential. The 
dependence of the KS potential on the semilocal current 
density of the system introduces a history-dependence 
not reflected by the ALDA or other adiabatic function- 
al of the charge or current density. The current density 
error Aj of a putative Kohn-Sham system that repro- 
duces the exact initial charge, but not current, density 
remains present in the KS potential thereafter via the lo- 
calized KS electric field which propagates alongside the 
wavepacket. 

The self-consistent method for calculating the exact 
Kohn-Sham electric field, while demonstrated for a one- 
dimensional system in Figure [2j is itself independent of 
dimensionality. While it is trivial to include such a field in 
the time-dependent Kohn-Sham equations in one dimen- 
sion, one could generalize to two- and three-dimensional 
systems with the inclusion of a KS vector potential, and 
an associated KS magnetic field, as well as a scalar poten- 
tial. Moreover, the method is, in principle, both universal 
and exact in the continuous limit, providing exact time- 
dependent Kohn-Sham potentials in the study of new, 
better functionals. 

The exact TDDFT Kohn-Sham potential for localized 
quasiparticle systems demonstrates the failure of local 
and adiabatic functionals such as the ALDA to describe 
them. A component of the exact potential has a func- 
tional dependence on the charge and current density that 
is highly nonlocal in both space and time, and is a cru- 
cial component of the complementary mechanisms that 
the exact potential employs. However, the inclusion of 
such effects in explicit scalar-potential functionals is in- 
trinsically problematic. An equivalent problem arises in 
the velocity gauge, with the presence of the nonadiabatic- 



ity of the exact Kohn-Sham vector potential functional 
|16j . The explicit inclusion of a Kohn-Sham electric field 
functional, however, may provide an additional ingredi- 
ent that allows nonlocal components of the potential to 
be implemented via calculation of a local or semilocal 
functional of the charge and current density. 

We thank Peter Bokes for fruitful discussions, and ac- 
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